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Abstract 

The successful anthropic prediction of the cosmological constant depends cruciahy on the as- 
sumption of a flat prior distribution. However, previous calculations in simplified landscape models 
showed that the prior distribution is staggered, suggesting a conflict with anthropic predictions. 
Here we analytically calculate the full distribution, including the prior and anthropic selection 
effects, in a toy landscape model with a realistic number of vacua, N ~ 10^°°. We show that it 
is possible for the fractal prior distribution we find to behave as an effectively flat distribution 
in a wide class of landscapes, depending on the regime of parameter space. Whether or not this 
possibility is realized depends on presently unknown details of the landscape. 
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I. INTRODUCTION 



The observed value of the cosmological constant A is about 120 orders of magnitude 
smaller than theoretically expected^ 



An ~ 10 



^120 



(1) 



One explanation for this observation assumes that A is an environmental parameter which 
has different values in different parts of the "multiverse" [l|, [s], 0, H, [g], 0, @] • The probability 
for a randomly picked observer to measure a given value of A can then be expressed as [3| 

Pobs(A) (X P(A)nobs(A), 



(2) 



where -P(A) is the prior distribution or volume fraction of regions with a given value of A 
and nobs (A) is the anthropic factor, which is proportional to the number of observers that 
will evolve per unit volume. If we assume that A is the only variable "constant", then 
the density of observers is roughly proportional to the fraction of matter clustered in large 
galaxies, nobs (A) oc /g(A). Using the Press-Schechter approximation ^] for /g(A) we can 
write [l0| 

, A ^ 

^obs(A) 



erfc 



A, 



(3) 



where we have normalized nobs to be 1 for A = 0. We have parameterized the anthropic 
suppression with a value Ac, which depends on such things as the amplitude of primordial 
fluctuations and the minimum size of galaxy that can contain observers. For the parameters 
used in Refs. 10, 11 , A^. is about 10 times the observed value of A, 



Ac ~ 



6 X 10 



-120 



(4) 



For the qualitative arguments and toy model of the present paper, the precise value of 
Ac will not matter, and we can use a Gaussian instead of erfc, to get 



"-obs(A) 



-(A/Ac)2''3 



(5) 



The prior distribution P(A) depends on the unknown details of the fundamental theory 
and on the dynamics of eternal inflation. However, it has been argued 12|, ll3| that it should 
be well approximated by a flat distribution. 



P(A) ^ const. 



(6) 



because the window where nobs (A) is substantially different from zero, is vastly less than the 
expected Planck scale range of variation of A. Any smooth function varying on some large 
characteristic scale will be nearly constant within a relatively tiny interval. Thus from Eq. 

(ED, 

Pobs(A) oc nobs (A)- (7) 

Indeed the observed value of A is reasonably typical of values drawn from the distribution 
of Eq. ([3]). This successful prediction for A depends on the assumption of a flat volume 



^ Here and below we use reduced Planck units, A/ijp = Mp/Sn = 1, where Mp is the Planck mass. 
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distribution ([6]). If, for example, one uses -P(A) oc A instead of ([6]), the 2cr prediction would 
be A/Aq < 500, giving no satisfactory explanation for why A is so small jl4|. 

A specific class of multiverse models is given by the landscape of string theory [l5, 16, 17 . 
In such models there are of order lO^'^'' different vacua with various cosmological constants 
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|20|. The dynamics of eternal inflation populates the multiverse with all possible 
vacua (or bubbles) by allowing for the nucleation of one vacuum within the other, according 
to transition rates which determine the probability of going from one vacuum to another. 

Given a specific string theory landscape, we would like to be able to predict the cosmo- 
logical constant that we should expect to observe according to Eq. ([2]). We will assume the 
prior probability P(A) is given by the relative bubble abundances of different vacua. Since 
an eternally inflating multiverse contains an infinite number of each type of vacuum allowed 
in the landscape, it is necessary to use some regularization procedure to compute the prior 
probability distribution. Many such regularization procedures, or probability measures, have 



23, 24, 25 



been proposed [2l|, [2 
13 and the references therein. 



Refs. \2 



Refs. 



30 



For a more complete and up to date account see Ref. 
Here we will use the pocket-based measure introduced in 
computed prior probabilities^ of different vacua in toy models 



15, l3] and did not find a smooth distribution of possible cosmological constants. Instead, 
for the specific models and parameters they studied, there were variations of many orders 
of magnitude in the prior probabilities of different vacua. However, to allow for numerical 



e p 

solution, Refs. [29|, 30| used models with a relatively small number of vacua and worked only 
in a first-order approximation. 

Here we will consider a toy model in which transition probabilities are computed as though 
all changes in A were by some fixed amount c. In this simple model, we can analytically 
study probability distributions for a realistic number of vacua, N ~ 10^°°. We find that 
when c is around 1, there is a smooth distribution of vacua in the anthropic range, and the 
anthropic prediction of Eq. ([7j) applies. But when c is smaller by a few orders of magnitude, 
the behavior is very different. In this case, the -P(A) factor is more important than ?T,obs(A) 
in Eq. ([2]). Thus we would expect to live in a region with large A, and so only a few galaxies. 
In such a multiverse, the anthropic procedure would not explain the observed small value of 
A. 

The plan of this paper is as follows: In section [TTl we will outline the method used to 
calculate bubble abundances. Additional details of the bubble abundance calculation are 
presented in the appendix. We will then define our model in section IIIIl and calculate the 
prior probability distribution. In section HVl we will investigate the behavior of Po6s(A). We 
end with a discussion in section |Vl 



II. BUBBLE ABUNDANCES 



We review here the procedure for calculating the volume fraction of vacua of a given kind, 
using the "pocket-based measure" formalism of Refs. 27, 29l . 



Vacua with A < are said to be terminal. There are no transitions out of them. Vacua 
with A > are recyclable. If j labels such a vacuum, it may be possible to nucleate bubbles 
of a new vacuum, say i, inside vacuum j. The transition rate Hij for this process is defined 
as the probability per unit time for an observer who is currently in vacuum j to find herself 



^ Strictly speaking bubble abundances were calculated. 
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in vacuum i. Using the logarithm of the scale factor as our time variable, 



K,,=r,,-H-\ (8) 

where Tij is the bubble nucleation rate per unit physical spacetime volume (same as \ij in 
0) and 

H, = (A,/3)V2 (9) 

is the expansion rate in vacuum j. 

Transition rates depend on the details of the landscape. However, if Aj < Aj, the rate of 
the transition upward from i to j is suppressed relative to the inverse, downward transition, 
by a factor which does not depend on the details of the process j3l[ |. 



-247r2 



1 1 



(10) 



Given the entire set of rates Kij, we can in principle compute the bubble abundance Pa for 
each vacuum a, following the methods of Refs. 27^2^- An exact calculation would require 



diagonalizing an N x N matrix. But as in Ref. [217 we can make the approximation that 
all upward transition rates are tiny compared to all downward transition rates from a given 
vacuum (see also the appendix). In that approximation, we can compute probabilities as 
follows. 

First, define the total down-tunneling rate for a vacuum j, 

D, = J2 ''v ■ (11) 

Ai<Aj 

Then define the dominant vacuum, referred to as vacuum *, as that recyclable vacuum 
whose Dj is the smallest. Since bubble nucleation rates are suppressed in low-energy vacua, 
we expect A^, to be fairly small, however we would not expect it to be so small as to be in 
the anthropic range. With physically sensible Hij it can be shown that this vacuum will have 
no downward transitions to vacua with positive A. To see that this is true, imagine that in 
some direction A* can jump downward to Aq, > 0. Now if we compare to we see that 
each term contributing to Da is less than the corresponding term (i.e., the transition rate 
in the same direction) in D^: because Aq, < A* and jump sizes in the same direction are the 
same. This implies Da < -D* which contradicts our definition of as the vacuum with the 
smallest sum of downward transition rates. 

Once we have identified the dominant vacuum, the probability for any vacuum a is given 
by (see Appendix) 

{Da-D,){D,-D,)---{D,-D,) ^ ^ 

where the sum is taken over all chains of intermediate vacua a,b, . . . , z that connect the 
vacuum a to the dominant vacuum. 



III. TOY MODEL 
A. Model 

We will consider a toy version of the Arkani-Hamed-Dimopolous-Kachru (ADK) model 



17 . We let there be 2J directions and so = 2^-^ vacua. We will choose J ~ 800, so that 
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TV ~ 10^°°. Each vacuum can be specified by a list of numbers {771, . . . , 772 j}, where rji — ±1, 
and the cosmological constant is 

^ = ^ + ^E^^^ (13) 

i 

The "toy" feature of this model is that all jumps have the same size, c. 

We will take the average cosmological constant A to be in the range (0, c). All vacua with 
J + coordinates and J — coordinates will have A = A. Each vacuum has 2J neighbors to 
which it can tunnel by bubble nucleation. Each nucleation event either increases or decreases 
the cosmological constant by c. 

The model as given above is of no use for anthropic reasoning. Vacua exist only with 
widely separated cosmological constants, A, A + c. . . , therefore we would not expect any in 
the anthropic range. So we will modify the model by artificially perturbing the A of each 
vacuum to produce a smooth number distribution. Vacua originally clustered at A will be 
spread out over the range from to c. This will cover the anthropic range of vacua with 
A > 0, and so if the vacua are dense enough we will find some anthropic vacua. We will not, 
however, take account of these perturbations in computing probabilities. 

We will only be interested in the vacua near A = 0, which are those at A before the 
perturbation procedure above. All these have exactly the same transition rates, and thus 
there is no single dominant vacuum. So, in addition to the "smearing" above, we will make 
a small perturbation to decrease the total tunneling rate of some specific vacuum A*, so that 
it is the dominant one and the procedures of the last section can be applied. 



B. Distribution of vacua 

The dominant vacuum, and all other vacua of interest have J + coordinates and J — 
coordinates. Thus the other vacua are reached from the dominant vacuum by taking an 

equal number of up jumps and down jumps. We will classify the vacua by a parameter n, 
the minimum number of up jumps required to reach a given vacuum from the dominant 
vacuum. We will call n the level of the vacuum. Thus a vacuum of level n differs from the 
dominant vacuum in 2n coordinates, n of which are + where the dominant vacuum had — , 
and another n vice versa. 

The total number of vacua of level n is thus 

We imagine these to be smeared over a range c, so their density is 

Pn^Njc=l/A^. (15) 

The likelihood that there is no vacuum in a range of size x is exp(— p„x), thus the median 
A of the lowest-A vacuum is 

A„ = (ln2)/p„ = c(ln2)/Ar„ (16) 

We will take a typical realization to be one whose lowest-A vacuum is at this median position. 
Above the lowest-A vacuum of level n, there are Nn — 1 more with higher A, with the typical 
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interval in A being A„. For a typical realization it is sufficient to take these vacua as evenly 
spaced, so that they are at 

Kn,j = A„ + (j - 1)A„ (17) 

where 1 < j < Nn- 



C. Probabilities 



We now use the formalism of Sec. [TTl to calculate the relative abundances of different 
vacua in our toy model. 

The relative abundance of each vacuum a is given by a sum over all chains that connect 
it to the dominant vacuum, Eq. (fT2l) . The minimum number of transitions in such a chain 
is 2n. Longer chains can be formed by jumping one way and then later the opposite way 
in the same direction. But these chains will have extra suppression factors because of the 
extra jumps, so the probability will be accurately given by including only minimum-length 
chains. Furthermore, the paths that maximize the bubble abundances are those that entail 
first making all the up jumps and then following with a sequence of down jumps. The reason 
is that the up-jump suppression factor, Eq. ffTOj) . is least when the starting A for the jump 
is highest. Thus it is best not to jump down until one has made all the necessary up-jumps. 
So we only consider the contribution of paths which consist of making all upward jumps 
first and then following with downward jumps. In this case, we can reorganize Eq. f|T2l) . 

E^aa l^ab ^^rs f^tu ^z* /'1Q^ 



D,Db-D, Ds-D, Dt 

The transition rates to the right of the factor in Eq. fll2l) are upward rates, and those 
to the left are downward rates. Hgt represents the first downward jump after having made n 
upward jumps from the dominant vacuum. 

Now we will approximate <^ Dj, since the transition rates are suppressed for low A 
vacua. Furthermore, in our single jump size model, all downward jumps from the same site 
have the same transition rate, so for Aj < Aj, 

^ (19) 



J. 



where Jj is the number of + coordinates in vacuum j. 
Using Eq. (fTOj), 



exp 



A, 



The product of all such suppression factors is just 



S = exp 



where A„ = nc is the maximum A reached. 



(20) 



(21) 



Factoring out the suppression factors, we find n terms given by Eq. (IT^ for the up jumps. 
The first one has Jj = J + 1, since there are J + coordinates in the dominant vacuum. The 
next has J + 2, and so on up to J + n. Thus the product is 



J!/(J + n)! 



(22) 
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The down-jumps are similar, except that the one at n is missing, giving 



J\/{J + n-l)\ 



(23) 



The up- and the down-jumps can be taken in any order, so there are (n!)^ equally weighted 
paths to reach the same vacuum. Thus the prior probability of a vacuum of level n is 



Pn cc {n\y 



_(J + n)! 



(J -|- n)Kst exp 



247r' 



1 1 

cn A, 



We set the first downward jump rate, which has no canceling denominator. 



Kst ~ exp 



a3/2 



exp 



(24) 



(25) 



which is the rate we would have for a Bousso-Polchinski(BP) type model 2^ with ^ AA 



IV. DISTRIBUTION FOR THE OBSERVED A 

Given the above, we are in a position to calculate the probability of observing each value 
Kn,j in a typical realization of our toy model. These are given by 



-fobs (-^n j) OC -Pj2?T.obs (An j) 

The chance that we live in a world of a given level n is then given by 

Pobs(?^) oc Pn ^ riobs(A„j) 



(26) 



(27) 



We will consider two cases. When A^ -C Ac, the sum can be approximated by an integral, 

3v/vrAc ?,^NnK 30rAcln2 

'0 



y'nobs(An,j) ~ -r- / (iAnobs(A) 



4A. 



4c 



4A, 



(28) 



Including Eqs. dH El [25]) , we find 

3v^Ac / J!^ 



Pobs('^) OC 



4c \{J - n)\{J + n)\ 



(J -|- n) exp 



247r2 eVSvr^ 

3/2 



cn 



cn 



(29) 



We have not included the term involving A^,, which is the same for all Pobs('^)- In Eq- 
( 129|) . Pobs(^) is a decreasing function of n. 

On the other hand, when A„ > A^ Eq. (1271) will be dominated by the first term, and we 
can write 

E^°bs(An,)^e-(^"/^=)''^ (30) 
j 

Including Eqs. ([211 ESD^ we find 



Pnhs(^) OC 



J\n\ 
{J + n)\ 



(J + n) exp 



37^" V a; 



cn 



cn 



(31) 
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In Eq. (13T!) . Pohs{n) increases with increasing n while n is small and the last term in the 
exponent is dominant, but it decreases when n is larger and the other terms are dominant. 
The division between regimes occurs when A„ ~ Ac, i.e.. 



J ^ ' 

n 



cln2 ~ Ac ~ 6 X 10"^^° (32) 



With c ~ 1, we find n ~ 34. The dependence on c is weak, with c ~ 10^^ corresponding to 
n ~ 33. For n in this range, changing n by one unit changes A„ by a factor of about 500. 
Thus there is at most one n with A„ ~ Ac. 

If we compare Eq. fl29l) and Eq. fl3T|) for the same n, we see that they differ by a factor 
of AcA^n/cexp (A„/Ac)^''^ ~ 1 if A„ ~ Ac, so there is no big jump due to switching regimes. 

Now let us start with n = 1 and increase n. Certainly with = 1, A,„ ^ Ac by a huge 
factor, we are in the regime of Eq. fl5T]) . and Pobs is infinitesimal. As we increase n, Pobs 
increases. Once n is significantly above 1, we can approximate the increase from one step 
to the next as 



Pobs(^ + l) rn 



Pohsin) VJ 



- exp 



247r2 ^ Qv^vr^ , /A, ^ 



cn 



2 2m5/2 VA 



(33) 



where we have ignored (A„+i/Ac)^'''^ as much less than (A„/Ac)^^'^. The ratio of the middle 
term in the exponent to the first term is {16y/n) ^ 0.05 for n ~ 33, so we will ignore 

the middle term. 

For sufficiently small n, the last term in the exponent dominates and Pobs (^+1) /-Pobs l'^) 3> 
1. There is only an infinitesimal probability that we will be in a vacuum of level n, because 
there are others that are much more probable. As we increase n, Pobs(^) will continue to 
increase. What happens next depends on the magnitude of c. 

A. Small c 

First suppose c is small, in particular that 

c < ^2/3j4/3 (34) 

for relevant values of n. For J = 800, n ~ 33, the right hand side is about 3 x 10""^. From 
Eq. m, 

247r2 

Cn2/3 J4/3 > 1 • (35) 

Now A„+i/A„ ^ (n/J)'^, and so successive values of An^ differ by a factor about n^l'^ j J^l'^ . 
Thus we can find a value of n such that 

247r2 /A„\^/^ 247r2 



We will now show that for this n, 

Pobs(n+l) <Pobs(n), (37) 
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so that we should find ourselves in a vacuum of at most level n. 

It is not clear from Eq. (136|) whether A„_|_i/Ac is more or less than 1, so we might need to 
use either Eq. ( I29l) or Eq. (!3T!) for Pobs('^ + !)• We will prove the claim using Eq. (l29l) . Since 
this gives a larger value than Eq. (PT]) . if Eq. (P7|) holds using Eq. (12^ . it will certainly hold 
using Eq. (pTi) . Thus we will take 



Pohsin) 



4A 



n+l 



Now from Eq. fl5Sl) . we find that 



j) exp 



247r2 ^ ^A^ 



247r2 ^ /jy/^ 



A, 



2/3' 



(38) 



(39) 



For J = 800, n ~ 33, the right hand side is about 70. Thus unless (A„/Ac)^/^ is extremely 
close to the upper bound in Eq. (l36l) . the exponential term in Eq. ( |38l) will be infinitesimal, 
and Eq. (EZD will follow. If we do have {K/ ^cf'^ ^ 247rV(cn2), then A„+i/Ae > 1. Then 
the prefactors in Eq. ( |38l) are at most about (n/ J)^, about 2 x 10"'^ for parameters of interest, 
and again Eq. fl37|) follows. 

Thus we can say with great confidence that we live in a universe with level n or lower. 
From Eq. fl36|) . we see immediately that we should observe A > A^, whereas in fact we observe 
Aq ~ O.lAc. If c is significantly smaller than the limit in Eq. flM|) . then we will see a very 
"non-anthropic" universe. We will be able to find A„ with (A„/Ac)^/^ > 247r^/(cn^/^ J"^/^), 
and thus 

(40) 



r^obs(A)<r^obs(An)<e-^^'^^/(^"^^^^*^^) 



will be tiny, meaning that only an infinitesimal fraction of matter has coalesced into galaxies. 
For example, with c = 10~^, we would find nobs (A) ~ ~ 0.05, in contrast to the 
observed value (in our approximation) ?7,obs(A) ~ 0.85. With c = 10"'^, we would find 

?^obs(A) 



< 



-30 



10 : a universe utterly unlike our own. 



B. Large c 

Now suppose instead that 

c > 247rVn^ (41) 

For J = 800, n ~ 33, the right hand side is about 0.2. Then we will reach the point where 
An ~ Ac before 247r^/ (cn^) (or the variation due to Hgt) is significant. In that case, we switch 
to the regime of Eq. (!29|l . where Pobs decreases only slowly with increasing n. In this regime, 

Pobs(n + l)^^ (J-nf ^.f^_!^Y (42) 



PoM {J + n + iy V J 

which is about 0.85 for parameters of interest. Thus we find that several values of n con- 
tribute nearly equally to the total probability. The first of these might be dominated by a 
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single A„, but the others will have a large number of closely spaced A. These vacua have 
similar nobs and identical prior probability, so we could easily be in any of them.^ 

Thus when c is large, we recover approximately the original anthropic predictions with 
a smooth prior -P(A). There might be an effect due to the discrete nature of the vacua 
associated with the smallest n, where Pohs{n) has its peak, but this effect is small because 
level n does not dominate the probability distribution. Instead the probability is divided 
across many different levels, while only level n has the above effect. 



V. DISCUSSION 

A key ingredient in the anthropic prediction of the cosmological constant is the assump- 
tion of a flat prior distribution. However, the first attempt to calculate this distribution 
for the Bousso-Polchinski and Arkani-Hamed-Dimopolous-Kachru landscape models 2^, 30l 
revealed a staggered distribution, suggesting a conflict with anthropic predictions. 

These calculations have been constrained by computational limitations and reveal only 
the probabilities of a handful of the most probable vacua^. In this paper we have gone 
beyond these first order perturbative results by studying a simple toy model which permits 
analytic calculation with a large, realistic number of vacua, ~ 10^°°. We have found 
an interesting fractal distribution for the prior -P(A). When including anthropic selection 
effects to determine Po6s(A), we find that agreement with observation depends on the only 
free parameter of the model, the jump size c. 

We have shown that when c ~ 1, anthropic reasoning does indeed solve the cosmological 
constant problem. Even though the prior distribution has a rich fractal structure, the states 
of interest have similar vacua sufficiently closely spaced to approximate the fiat distribution 
well enough to give the usual anthropic results. 

Bousso and Yang S^l discuss the probability distribution resulting from the pocket- 



based measure of Garriga et al [27j. They claim that it cannot solve the cosmological 
constant problem, because the Shannon entropy, S = — ^plnp, computed from the prior 
probabilities will never obey exp(S') ^ 10^^°. We feel that a better test is to compute the 
entropy using the same formula with the probabilities taking into account nobs, and then to 
demand that 5* 3> 1, so the effective number of places in the landscape in which we might 
find ourselves is large. This condition is clearly obeyed in the case where c ~ 1. 

On the other hand, if c is small, of order 10~^ for J = 800, then the agreement with 
observation breaks down. In this case, we should expect to find ourselves in a universe with 
quite a large cosmological constant. Even though the anthropic factor strongly disfavors 
such universes, their volume fraction is so much higher that in the overall probability they 
are greatly preferred. 

Of course these results apply directly to our toy model with nearly identical jumps, and 
there is no reason that the real landscape of string theory should have this property. What 



3 



4 



The identical prior probabilities are a toy feature of the model. But even if we were to distribute these 
probabilities between Pn and Pn+i, we would still find many vacua with similar probabilities. If level n 
has 1 vacuum that is not significantly suppressed by nobs, then level n + 1 will have Nn+i/^n such vacua. 
These vacua will be distributed in a range of probabilities with Pn/Pn+i ^ N„+i/Nn, so it is not possible 
to have a single one of them strongly dominant. 

It is extremely unlikely that any of these vacua should lie in the anthropic range. 
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would happen in a more realistic theory, which would have many different-sized jumps in 
A? 

Consider first what happens when we rescale our theory by changing all jump sizes by 
a constant factor. Suppose, for example, that all different jump sizes are proportional to a 
parameter c.^ The transition rates between vacua all have terms proportional to 1/c in the 
exponent. Thus when c is small, the transition rates and thus the probabilities of different 
vacua are very sensitive to the details of the transition rates. When c is large, all rates are 
larger and less variable. Thus we conjecture that, as a general rule, landscapes with large 
jumps are more likely to give the standard anthropic results, while those with small jumps 
are likely to predict universes unlike ours. 

Regardless of the details of the theory, there is always a factor proportional to 
exp(247r^/Amax)! where A^ax is the maximum cosmological constant reached in the series of 
vacua between the dominant vacuum and ours. If it is necessary to jump up to a high value 
of Ajnax in order to have a large enough number of possibilities to expect any vacuum near 
the anthropic range, then the exact value of A^ax will not be so important. If, however, 
enough vacua can be reached with very small Amax? then this factor will depend strongly on 
their exact values of Amax- 

If this effect is the dominant one (which is not clear in a realistic model), then what 
matters is not the average size of the jumps, but the number of small jumps. If there 
are enough small jumps to produce some near-anthropic vacua, then these vacua will be 
preferred, because of their small Amax, and there may be disagreement with observation. 
The existence of large jumps in addition is then not important. Work is underway to study 
models with different jump sizes. 



APPENDIX A: BUBBLE ABUNDANCES BY PERTURBATION THEORY 



As shown in Ref. [27|, the calculation of bubble abundances pj reduces to finding the 
smallest eigenvalue q and the corresponding eigenvector s for a huge N x N recycling tran- 
sition matrix R. Bubble abundances are given by 

Pj ^ ] H^KjaSg ~ ^ ] t^jaSg- (^1) 

a a 

where the summation is over all recyclable vacua which can directly tunnel to j. Ha is the 
Hubble expansion rate in vacuum a and we take H'^ ^ 1 because q is an exponentially small 
number. 



In a realistic model, we expect to be very large. In the numerical example of Ref. 29 1 

19, 



N ~ 10^, while for a realistic string theory landscape we expect ~ 10^°'' (igI . 
20[ |. Solving for the dominant eigenvector for such huge matrices is numerically impossible. 



However, in Ref. [29|, |33| the eigenvalue problem was solved via perturbation theory, with 
the upward transition rates (see Eq. (fTOj) ) playing the role of small expansion parameters. 
Here we extend this procedure to all orders of perturbation theory. 

We represent our transition matrix as a sum of an unperturbed matrix and a small 
correction, 

R = R(°) + R('), (A2) 



^ For example, consider a model with |AAi| = ci where 1 < i < 2 J. 2 J is the number of directions and c is 
still some overall scale. 



11 



where R^^'' contains all the downward transition rates and R*^-*^^ contains all the upward 
transition rates. We will solve for the zero'th order dominant eigensystem {g*^°\s'^^^} from 
R^"-* and then include contributions from R^-""-* to all orders of perturbation theory. 
If the vacua are arranged in the order of increasing A, so that 

Ai < A2 < ••• < Aiv, (A3) 

then R(°) is an upper triangular matrix. Its eigenvalues are simply equal to its diagonal 
elements, 

j<Ol 

Hence, the magnitude of the smallest zeroth-order eigenvalue is 

=D, = min{D„}. (A5) 

Downward transitions from * will bring us to the negative-A territory of terminal vacua 
29| . Terminal vacua do not belong in the matrix R; hence, i?/3* = for (3 ^ and we see 



that the zeroth order eigenvector has a single nonzero component, 

= (A6) 

Thus, in fact, we could have included only the diagonal elements of R'-^^ and still obtained 
the correct q^^^ and s*^°\ 

Now we would like to include the effect of the upward transition rates in the lower trian- 
gular matrix R*^°) , to any given order of perturbation theory. To organize the calculation, 
we note that the eigenvalues of an upper triangular matrix are just the diagonal elements, 
and the eigenvectors are given exactly by the perturbation series, which terminates after N 
terms. Thus we will take just the diagonal elements of R'-"^ as our unperturbed matrix, and 
consider the rest of R*^"^ and R'--'^^ as a perturbation. By summing all terms in the pertur- 
bation series for this perturbation that involve n elements of R*^"'"^ we find the perturbation 
term of order n in the small upward jump rates. 

The procedure is the same as that involved in finding an eigenstate of the Schrodinger 
equation in nth-order perturbation theory, working in a basis of wavefunctions which diag- 
onalize the unperturbed equation. See for example Eq. (9.1.16) of Ref. [3J]. Including all 
orders of perturbation theory, we find the result 

N N 

where there can be any number of terms in the sum and the vacuum * is not summed over. 
Combining Eqs. (jMD and ([ATD gives Eq. ffT2D . 
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